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A null ray approaching a distant astronomical source appears to slow down, while a 
massive particle speeds up in accordance with Newtonian gravitation. The integration 
of these apparently incompatible aspects of motion in general relativity is due to the 
existence of a critical speed. Dynamics of particles moving faster than the critical speed 
could then be contrary to Newtonian expectations. Working within the framework of 
gravitoelectromagnetism, the implications of the existence of a critical speed are explored. 
The results are expected to be significant for high energy astrophysics. 
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1. Introduction 

Gravitoelectromagnetism ("GEM") originates from the similarity between 
Coulomb's law of electricity and Newton's law of universal gravitation. In the course 
of the development of his dynamical theory of the electromagnetic field, Maxwell 
considered a similar theory for the gravitational field 1 . He noted that gravitational 
charges are all of the same kind and lead to attraction instead of repulsion. On the 
basis of certain energy arguments, he then concluded that one could not arrive at a 
fundamental theory of gravitation in this way 1 . Several years later, difficulties with 
the excess perihelion motion of Mercury led to the phenomenological introduction 
of a gravitomagentic force due to mass current in the Sun 2 ' 3 . However, in 1915 
Einstein's general relativity provided a beautiful explanation of Mercury's motion 4 . 
This involved a relativistic correction to the Newtonian gravitoelectric potential of 
the Sun. Soon afterwards, the gravitational influence of the rotation of the Sun on 
planetary orbits was determined within general relativity by de Sitter 5 . This equa- 
torial gravitomagnetic precession turned out to be much smaller and in the opposite 
sense as the Einstein gravitoelectric precession 5 . 

A field theory that successfully combines Newtonian gravitation with Lorentz 
invariance should necessarily contain a gravitomagnetic field 6 . The first general in- 
vestigation of the gravitomagnetic field within the framework of general relativity 
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is due to Thiriing 7 . Moreover, the general rate of precession of the orbit of a test 
particle in the field of a rotating mass was determined by Lense and Thirring 8 . An 
English translation of the Thirring-Lense papers together with a critical commen- 
tary on their results has been provided in Ref. 9. Extensive reviews and discussions 
of the various aspects of GEM are contained in Refs. 10 -12. On the observational 
side, there has been progress, especially since NASA's GP-B was launched on April 
20, 2004. 13 " 15 

The GEM framework will be employed in this paper in order to illustrate the ex- 
istence of a critical speed in the dynamics of a test particle moving in a given GEM 
field. There are two principal approaches to GEM. The linear perturbation approach 
to GEM is briefly presented in Section 2 and the general equation of motion of a 
test particle is investigated in Section 3. It is shown that this equation of motion 
contains a critical speed given by v c = c/y/3 in the case of one-dimensional motion. 
Section 4 is devoted to a more basic invariant treatment of the critical speed; for this 
purpose, the spacetime curvature approach to GEM is employed using the quasi- 
inertial Fermi coordinates. In this invariant analysis of relative motion, the critical 
speed turns out to be V c — cj\[2. For relative ( "ultrarelativistic" ) motion that is 
faster than this critical speed, the results of recent investigations of tidal accelera- 
tion/deceleration phenomena in the field of a neutron star or a black hole are briefly 
summarized. The case of ultrarelativistic fluid flow is very similar and the relevant 
equations are derived in Appendix A. Finally, Section 5 contains general remarks 
regarding the gravitational acceleration/deceleration of ultrarelativistic particles. 



2. Motion in a GEM Field 

Gravitoelectromagnetism provides a useful method for the description of the gravi- 
tational field generated by a slowly rotating "nonrelativistic" astronomical source in 
the linear approximation of general relativity 4 . The spacetime metric can be writ- 
ten as g^ v = r]^ + hfjiv, where rj^ u is the Minkowski metric tensor with signature 
+2 and is a first-order perturbation. The gravitational potentials are in gen- 
eral gauge-dependent, i.e. — > h^ v + e^„ + e„ iM , where e M is a vector associated 
with the choice of the background inertial coordinates = [ct, x). In terms of the 
trace-reversed potentials = h^ u — ^rj^h, where h = tr(/i^), the gravitational 
field equations can be written as 

nl 16TrG T m 
L-lh^ = W 

after the ( "Lorentz" ) gauge condition W v v = has been imposed. 
The general retarded solution of (1) involves the particular solution 



h f 

^ ~ C 4 J 



V(ci-|x-x'|,x') , 3 



x — X' 



dV (2) 



plus a solution of the homogeneous wave equation that we simply ignore in the 
present analysis. For the sources under consideration, hoo = 4<I>/c 2 , hoi = —2Ai/c 2 
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and hij = 0(c~ 4 ). Here x) is the gravitoelectric potential and A(t, x) is the 
gravitomagnetic vector potential. All terms of 0(c~ 4 ) will be neglected in the metric 
tensor. Thus the important quantities of interest in equation (2) are T 00 = pc 2 and 
T° l = cj l , where j 7 * = (cp,j) is the mass-energy current. The conservation of the 
mass-energy of the source is assured through the imposition of the gauge condition, 
i.e. j M =0 follows from 

£ + v.(iA)-a 

The spacetime metric is thus given by 

-ds 2 = -c 2 (\ - 2^ dt 2 - ~ c {A ■ dx)dt +(^ + 2 J) Sijdx'dx*. (4) 

It is possible to define the gravitoelectric and gravitomagnetic fields in close 
analogy with electrodynamics 

e =- v *4!G a )' b = vxa - (5) 

It follows from equations (3) and (5) that 

VxE = -I|fV), V.fVUo, (6) 



cdt\2 J' \2 
while the gravitational field equations (1) reduce to 

V-E = 4^, Vx( J Bj=--E + _ J . (7) 

These are the Maxwell equations for the GEM field; they are based on a certain 
convention that is explained in the following paragraph. 

To preserve the electromagnetic analogy as much as possible, it would be conve- 
nient to be able to employ in the GEM case the standard results of electrodynamics 
using a special convention. To this end, we assume that the source has gravitoelec- 
tric charge Qe = GM and gravitomagnetic charge Qb = 2GM, where M is the 
total mass of the source. It follows that the gravitomagnetic dipole moment of the 
source is G3 /c, where J is the total angular momentum of the source. Thus far from 
the source (r = |x| > GM/c 2 ) 

^ GM GJxx 

' A ~7^^' ( 8 ) 

which are consistent with equations (2) and (3). Moreover, a test particle of inertial 
mass m has gravitoelectric charge qE — —m and gravitomagnetic charge qs = —2m 
in this convention. We note that the signs of (qE>QB) are opposite to those of 
(Qe,Qb) to preserve the attractive nature of gravity; furthermore, the ratio of 
gravitomagnetic charge to the gravitoelectric charge is always 2, since the linear 
approximation of general relativity involves a spin-2 field. This is consistent with 
the fact that the ratio of the magnetic charge to the electric charge of a particle is 
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unity in Maxwell's spin-1 electrodynamics 16 . The magnetic charge employed here 
should be distinguished from the magnetic monopole strength, which is strictly zero 
throughout this paper. 

The geodesic equation of motion of a test particle in the stationary GEM field 
of a source with d<fr/dt = and dA/dt = may be written in a form analogous to 
the Lorentz force law, 

mT = -mE-2m-xB, (9) 

at c 

when velocity-dependent terms of order higher than v/c are neglected. To go be- 
yond the GEM analogy, let us explore the significance of the terms that have been 
neglected in equation (9). The geodesic equation for a particle with proper time r 
and four- velocity — dx^/dr — 7(1, v/c) is given by 

~dV 

Once the linear GEM field is given, it is possible to employ equation (10) for the 
motion of particles of any speed in this field; for instance, one may consider the 
force-free motion of cosmic rays in the gravitational field of the Earth. For the sake 
of simplicity, we choose units such that c = 1 in what follows. For the GEM field, 
the Christoffcl symbols are 

rg„ = -*,„, r?. = 2A (<J) + <Sy$, , (11) 

r 00 = - 2A ifi , r Qj = 5y$ i0 + e ijk B k , (12) 
v )k = fc + <Sik*j - fyfe* i- (13) 

The physical consequences of the geodesic equation (10) with the connection given 
by equations (11)-(13) are explored in the following section. 

3. Critical Speed 

It is straightforward to show that the components of equation (10) can be expressed 
as 



ld7 - = ( l_ r 2 

dv 



= (l-v 2 )$ + 2v i [$ i -A (itj) v*], (14) 



7 dt 

(1 + v 2 )<S>. % - 2(v x B)i + 2A lfl - «*(3 - « 2 )$, 



dt 

+ 2v i vi[A {m v k -2$ J ]. (15) 

Here we have used the fact that du^/dr = jdu^/dt and then separated the equations 
for 7 and v. The worldline of the test particle is timclike; therefore, it follows from 
u^u^ = — 1 that 

^ = l-i) 2 -2(l|» 2 )$ + 4vA. (16) 
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Consider a stationary axisymmetric source such that for r » GM/c 2 the po- 
tentials are given by equation (8). In this case, equation (15) may be expressed 
as 

dv GM,. o, A . , n 2G, 2t 
_ = -_[(! + v 2 ) x _ 4 (x • v)v] - — [r 2 J x v 

-3(J-x)L + 3(J-L)(x- v)v], (17) 

where L = x x v. For one-dimensional motion along the rotation axis of the source, 
equation (17) reduces to 

dv GMx „ 2 , 

This equation contains a critical speed v c = l/\/3; that is, for motion with v < v c , 
we have the standard attractive force of gravity familiar from Newtonian physics, 
while for v — v Cl the particle experiences no force and for v > v c the gravitational 
attraction turns to repulsion. These results are valid in the linear approximation 
for the gravitational field under consideration here. It is interesting to note that if 
we use instead the standard post-Newtonian approximation scheme for the field of 
the source, the factor (l — 3v 2 ) in equation (18) becomes (1 — 3v 2 — 4$), where 
$ = GM/r < l. 17 

Equation (18) has the exact solution (v — dr/dt) 

v 2 =v 2 c -(v 2 c -vl)e^; (19) 

however, in the linear approximation, e~ 6 * ~ 1 — 6$, and equation (19) is valid 
only in the form 

v 2 =vl + 2 -^(l-3vl), (20) 

so that the speed of an infalling test particle increases, remains constant or decreases 
depending on whether Voo is less than, equal to or more than v c — 1/ \/3, respectively. 
In the nonrelativistic limit, 1 — 'iv 2 ^ — > 1 and equation (20) reduces to the standard 
result of Newtonian gravitation. 

It is crucial to recognize here that highly relativistic outflows in the form of 
astrophysical jets are expected to occur along the rotation axis of the central massive 
source. Thus an equation such as (18), though valid for motion along the rotation 
axis, could nevertheless play a rather significant role in the gravitational physics of 
phenomena associated with jets. 

In connection with jets, it is useful to consider radial motion in the exterior Kerr 
spacetime. The Kerr metric is given by 

-ds 2 = -dt 2 + S ^dp 2 + dd 2 ^j + (f + a 2 ) sin 2 ddcf 2 

+ 2GM^(dt-asm 2 $d(t>) 2 , (21) 
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where £ = p 2 + a 2 cos 2 •& and A = p 2 — 2GMp + a 2 in Boycr-Lindquist coordinates 
(t, p, tf, <p). Here M is the mass and a = J/M is the specific angular momentum of 
the source. The geodesic equations of motion for a test particle moving along the 
rotation axis reduce to 



dj _ . p 2 + a 2 dp _ / 2 2GMp 

dr 1 p 2 -2GMp + a 2 ' dr Y p 2 + a 2 ' { ' 

where 7 is an integration constant. When 7 > 1, it is the Lorentz factor of the 
test particle as measured by the static inertial observers at spatial infinity i.e. 
7~ 2 = 1 — v 2 ^. It follows from system (22) that 

dpV 2 ( i~G 2 M 2 \ 2GM ( , 4GM\ n ( \\ 



dt) 00 v p 2 ) p V 00 p ) Vp 3 , 

To compare with our GEM results, this equation must be expressed in terms of the 
isotropic radial coordinate r, 



i + ^rj • (24) 

Thus equation (20) can be recovered, in terms of either p or r, to linear order in 
the gravitoelectric potential. 

Shapiro's radar echo delay experiments have demonstrated that photons in effect 
slow down in the gravitational field of a mass M . In Newtonian gravity, however, 
massive test particles speed up as they approach a gravitational source. It turns out 
that this is also the case in general relativity if the speed of the particle at infinity is 
below a critical speed. Otherwise, as Vao approaches the speed of light the particle 
has lightlike behavior. That a particle with an initial speed Voo above a critical 
speed slows down in the Schwarzschild field as the particle radially approaches the 
source was first demonstrated in Refs. 18 and 19. The critical speed v c = l/\/3 
was discussed by Carmeli 18 , who derived equation (20) for radial motion in the 
exterior Schwarzschild field. The propulsion aspects of this result have recently 
received attention. 20 Using a more invariant approach, Jaffe and Shapiro 19 discussed 
a transition velocity, vt — I/a/2, in connection with the radial motion of particles in 
the Schwarzschild field. They argued that the coordinate speed should be replaced 
with dl/dt, where dl 2 — j i jdx l dx : > and 7^ = gij — gmgoj/goo- This results, working 
to linear order in the Schwarzschild spacctime, in an equation of the same form as 
equation (20) except that the factor (1 — Sv 2 ^) is replaced by (1 — 2^). 

It is clear that by employing physically admissible coordinate systems one can 
obtain different numerical values for the critical speed. The critical speeds 1/V3 and 
I/1/2 have also been discussed in Ref. 21. Nevertheless, it is important to remark 
here that with respect to the basic class of static (generally noninertial) observers 
in the exterior Schwarzschild spacetime, the local speed of an infalling particle 
monotonically increases toward unity, while the local speed of light is always equal 
to unity. 
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In the standard post-Newtonian approximation scheme, one is limited to the 
slow-motion weak-field regime and therefore the coordinate quantities in the equa- 
tions of motion arc assumed to be physically significant as they are observable 
in the standard interpretation. Here, however, we have been using the linear post- 
Newtonian scheme for the determination of the field generated by the central source, 
while the test particles can travel at any speed. Thus to study further the predic- 
tions of general relativity in this case, we must construct appropriate invariants that 
would represent actually measurable quantities. To this end, we find it convenient 
to study relative motion in the quasi-inertial Fermi coordinate system 22 established 
about the motion of a reference observer. This manifestly invariant approach is 
described in the next section. 

4. Fermi Coordinates 

Consider an observer O following a geodesic in the gravitational field of an astro- 
nomical source. Let ^ be the observer's local orthonormal tetrad frame that is 

(<*) 

parallel propagated along the worldlinc of O. That is, ^ ^ = dx^/dr is the ob- 
server's four- velocity vector as well as its local temporal axis and i = 1, 2, 3 
are the spatial unit gyro directions that form the local spatial frame of the observer. 
Here t is the proper time of O. To describe physical phenomena relative to O, it 
is natural to construct a Fermi coordinate system along its worldline. This is a 
geodesic coordinate system based on the observer's local frame. Indeed, an event P 
with Fermi coordinates = (T, X) can be orthogonally connected to the worldlinc 
of the reference observer at Pq by a unique spacelike geodesic of proper length a. Let 
t be the proper time along the reference trajectory at Po and rf = (dx^ / da)o be the 
unit tangent vector to the spacelike geodesic at P : then, T = t and Xi = arj^X^^. 
Thus, O has Fermi coordinates (r, 0), so that the reference observer always occupies 
the spatial origin of the Fermi coordinate system. The spacetime metric in Fermi 
coordinates is given by 

F 3oo = - 1 - F RoiOj (T)X l X^ + . . . , (25) 
F g 0i = ~\ F R ajlk (T)X'X k + ..., (26) 

F gij = kj - \ F R tk] i (T)X k X l + . . . , (27) 

where 

F Raf3is(T) = R^vpa^ ( a )A" (/3) A p ((5) (28) 

is the projection of the Riemann tensor on the local tetrad frame of observer O. 
The Fermi coordinates are admissible in a cylindrical region with |X| < 1Z along 
the reference worldline such that 1Z(T) is a certain minimum radius of curvature of 
spacetime. 
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The equation of motion of a free test particle in the Fermi coordinate system 
can be expressed as 

" + F Kd Uau0 = 0. (29) 



ds af3 

where — dX^/ds — T(l, V) is the particle's four-velocity. The geodesic motion 
of the particle is always timelike; therefore, 

T- 2 = - F g o - 2 F g M V i - '*//,, I ' I •*' > 0. (30) 

It follows from equation (29) that 

ldT _ F dX^dX^ 

TdT af) dT dT ' ( ' 

The equation of relative motion (32) can be expressed to linear order in distance 
away from the reference observer O by using the terms given explicitly in the metric 
tensor (25)-(27). The result is 

^ + F R 0i0j Xi + 2 F R lkl0 V k Xi 

+ ^(3 F R okj0 V i V k + F R lkol V k V l + F R k3iV l V k V l )Xi = 0, (33) 
while the modified Lorentz factor is given by 

^ = 1 - V 2 + F R MQj X l X> + i F R 0jlk X^V l X k 

+ \ F R ikj iV i X k V j X l . (34) 

Neglecting the velocity-dependent terms in equation (33), one recovers the Jacobi 
equation. It is important to note that the treatment of test particle motion can be 
generalized to the flow of a perfect fluid as in Appendix A. 

To interpret the generalized Jacobi equation (33), let us first proceed by analogy 
with the linear perturbation theory of Section 2 and note that the metric in the 
Fermi frame reduces to the Minkowski metric along the reference trajectory; hence, 
one may write F goo = — 1 + 2$ and F got = —2A i: where 

$ = -\ F R 0t03 X*X\ A i = ^ F R 0jlk X^X k . (35) 

Moreover, the corresponding GEM fields are obtained from equation (5), so that to 
lowest order 

Ei = F R i0jX j , Bi = —-eijk F RjkoiX l . (36) 
It follows that equation (33) can be written as 

m ^- = qE E + qB VxB + 0(V 2 ), (37) 
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where qE = —m and qs — —2m as before. Though this analog of the Lorentz force 
law is extremely interesting and leads to the physical interpretation of the Bel and 
Bel- Robinson tensors 23 , it is important to recognize that the higher-order terms in 
equation (37) become rather significant for ultrarelativistic motion. To see this, let 
us imagine a circumstance where one-dimensional motion — along the Z direction, 
say — is possible. Then it follows from the symmetries of the Riemann tensor that 
equation (33) reduces to 

g + fc(T)(l-2Z 2 )Z = 0, (38) 

where Z = dZ/dT and k(T) = F Rtztz- The generalized Jacobi equation (38) has 
solutions for Z = ±l/\/2 such that the relative motion within the Fermi frame is 
uniform. Below the critical speed V c = l/V%, the relative motion is similar to what 
is expected from relativistic tides in accordance with the Jacobi equation. The sit- 
uation changes drastically, however, for ultrarelativistic relative motion above the 
critical speed V c . If k(T) < 0, then an initially ultrarelativistic particle decelerates 
with respect to the reference particle and asymptotically approaches uniform mo- 
tion with critical speed V c . If k(T) > 0, then an initially ultrarelativistic particle 
accelerates with respect to the reference particle. 24 These results for the motion of 
free particles can be extended to fluid flow (see Appendix A). 

Ultrarelativistic particles are expected to be produced in abundance in the cen- 
tral engines of active galactic nuclei as well as supernova explosions and X-ray 
sources. It is therefore natural to expect that the tidal acceleration/deceleration 
of ultrarelativistic particles could have significant implications for the physics of 
astrophysical jets and high-energy cosmic rays. 

Assuming that the gravitational field of the central source may be described by 
the Kerr field, consider a Fermi coordinate system established along reference escape 
trajectories on the axis of rotation given by equation (22) with 7 > 1. Such particles 
with 7 near unity provide an ambient medium surrounding the central source. High- 
energy particles originating near the poles of the central collapsed configuration and 
moving outward along the rotation axis faster than the critical speed l/y/2 w 0.7 
relative to the ambient medium experience significant tidal deceleration 25 , while 
such ultrarelativistic particles moving outward normal to the rotation axis expe- 
rience significant tidal acceleration relative to the ambient medium. In a series of 
recent papers, these phenomena have been investigated in some detail. 24-27 It turns 
out that tidal deceleration occurs for outflow velocities within a cone of half angle 
6 w tan -1 y/2 w 55° around the rotation axis. 26 ' 27 On the other hand, tidal accel- 
eration occurs for outflow velocities outside the critical velocity cone. 26 ' 27 

A complete analysis must of course include plasma and radiation effects as 
well. Nevertheless, preliminary results appear to be consistent with observations 
of microquasars and Chandra studies of certain neutron star X-ray sources in our 
Galaxy. 25 " 27 
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5. Discussion 

Using the main approaches to GEM , we have demonstrated the existence of a crit- 
ical speed that could have fundamental consequences for the gravitational physics 
of ultrarelativistic flows. The physical origin of the critical speed can be traced back 
to a "conflict" between the intuitively expected behavior of particles based on New- 
tonian gravitation and the behavior of light in a gravitational field; indeed, both of 
these aspects are properly integrated within general relativity. 

Imagine a light ray propagating from event Pi : (ti,xi) to event Pi : ^2,^2) 
in a global inertial system of reference; hence, f 2 — ii = |x 2 — xi|. If the spacetime 
is now slightly perturbed by the presence of a GEM field as in Section 2, then the 
arrival of the signal at P 2 is accompanied mainly by a Shapiro gravitoelectric delay 
given by 

A GE = 2 $ dl, (39) 
JPi 

where d£ is the element of length along the path of the unperturbed ray. 28 This 
is intimately related to the fact that the propagation of electromagnetic waves in 
the gravitational medium takes place with an effective index of refraction given 
by 29 n s=s 1 + 2$ for $ << 1. On the other hand, the kinetic energy of a particle 
increases as it gets closer to a gravitational source according to the Newtonian theory 
of gravitation. The conflict is resolved in general relativity through the existence 
of a critical speed such that "ultrarelativistic" particle motion with speed above 
this critical speed has significant features that violate our nonrelativistic intuition 
based on Newtonian gravitation. In the linear post-Newtonian approach, the critical 
coordinate speed is l/\/3; however, for relative motion in Fermi coordinates, the 
invariant critical speed is 1/V2. The latter case has been the subject of recent 
studies due to its direct observational significance. 

The most significant consequence of the existence of the critical speed is a pos- 
sible gravitational mechanism for the acceleration/deceleration of ultrarelativistic 
particles. This mechanism can be elucidated by comparing and contrasting the grav- 
itational and electromagnetic interactions. The motion of a charged particle in an 
external electromagnetic field in Minkowski spacetime is given by 

in accordance with the Lorentz force law. This equation may be written as 

Tt = ^^-^ [£ V(V ' £) + V X B] ' (41) 
which indicates that the closer v gets to unity, the more difficult it will be to change 
the velocity of the particle. For a finite external electromagnetic field, the motion 
of a charged particle with m 7^ is effectively uniform for v approaching unity. 
Thus the critical speed that appears in the electromagentic case is the fundamental 
speed, namely, unity. 
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The gravitational analog of equation (40) in the quasi-inertial Fermi coordinate 
system is 

^(rn = -r( F n Q + 2 F r 0j v= + F r jk vw k ). (42) 

while the analog of equation (41) is given by equation (32). The fact that the 
right-hand side of equation (42) is proportional to the modified Lorentz factor T 
is responsible for the bending of ultrarelativistic particle orbits by a gravitational 
field — indeed, this deflection approaches that of light in the appropriate limit. Sim- 
ilarly, equation (32) implies that it is in general possible to change the velocity 
of ultrarelativistic particles considerably in a gravitational field, since the analog 
of vl — v 2 in equation (41) involving the fundamental speed is missing in equa- 
tion (32); indeed, there is a critical speed of l/y/2, but it appears in a different 
way. An even more striking illustration of this circumstance can be given for one- 
dimensional motion. 

Consider a situation where £ and B are parallel and the charged particle moves 
along the field lines. Then equation (41) reduces to 

f.id-^e, ,43) 

which may be compared and contrasted with equations (18) and (38). While the 
ultimate ( "critical" ) speed in equation (43) severely restricts the magnitude of ac- 
celeration of a particle with v — > 1, this is not the case for gravitational motion in 
equations (18) and (38), which exhibit instead critical speeds of l/\/3 and l/y/2, 
respectively. 

The appearance of the critical factor (1 — 2Z 2 ) in equation (38) is due to the 
fact that in the quasi-inertial Fermi coordinates, one finds that to lowest order 

F V 00 = ^ = -1(^00),*- (44) 

Indeed, the first term in the critical factor originates with f Fq , while the second 
term originates with F Tg i . The factor 2 in the latter term is due to the fact that 
F r[] i occurs twice, since F T < Q i = F r" ; hence, the tensorial character of general 
relativity is responsible for the appearance of the factor of 2 in (1 — 2Z 2 ). Thus the 
circumstance that the critical speed is unity in electrodynamics while it is 1/V2 
in gravitation is related to the spin-1 nature of the electromagnetic field and the 
spin-2 nature of the linearized general relativity in the Fermi system, respectively. 

Recent investigations of the novel gravitational mechanism for the accelera- 
tion/deceleration of ultrarelativistic particles indicate that this mechanism may 
play a fundamental role in high-energy astrophysics. 24-27 Furthermore, the tidal 
motion of charged particles has also been recently examined in an interesting astro- 
physical context in Ref. 30, where the important role of the critical speed has been 
independently confirmed. 
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Appendix A. Gravitational Dynamics of Relativistic Flows 

Imagine a perfect fluid with energy density p and pressure p in the gravitational 
field considered in Section 4. The energy-momentum tensor of the fluid is given by 

T" v = (p + p)uV+pg^, (A.l) 

where w M is the unit four- velocity vector of the fluid. The purpose of this Appendix is 
to express the dynamical laws T^ v . v = in the Fermi coordinate system established 
about the worldlinc of the geodesic observer O. 

Let us first note that in a global inertial frame of reference, the dynamical laws 
may be expressed as 

where p — ~/ 2 (p + p) and 7 = (1 — v 2 )^ 1 / 2 is the Lorentz factor of the fluid. 
Equation (A. 2) follows directly from the p = component of T^ v . v — 0, while 
the p = i component, when combined with equation (A. 2) implies equation (A. 3). 
Equations (A. 2) and (A. 3) are the relativistic continuity and Euler equations, re- 
spectively, and assume their familiar form when the pressure has no explicit depen- 
dence upon time. 

In the quasi-incrtial Fermi coordinates, we expect that the dynamical laws would 
have the same form as equations (A. 2) and (A. 3) but modified by the presence of 
curvature terms. In fact, let the four-velocity be T(l, V) as before, where T > is 
given by equation (30). Then with p = T 2 (p + p), the modified equations of motion 
become 

|£ + V • (pV) = -p, Q V" - pK, (A.4) 



dT +{ a ' 3 afjV > dT dT 



V" + (P,a F 9 0a )V\ (A.5) 



where 



To lowest order, K is given by 

K = 1(8 F R 0i0j - F R l3 )V*XJ + 1 F R 0ijk V i V^X k . (A.7) 

It follows from these results that equations (A. 2) and (A. 3) are recovered along 
the reference trajectory at (T, 0), as required. Moreover, in the absence of pressure 
equation (A.5) is equivalent to equation (32) with attendant consequences regarding 
the existence of the critical speed V c = l/\/2. A complete analysis of equations (A.4) 
and (A.5) is beyond the scope of this work. 
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